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Renormalized QCD-inspired model for the pion and mesons 
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We apply the subtraction method to an effective QCD-inspired model, which includes the Coulomb plus a zero- 
range hyperfine interactions, to define a renormalized Hamiltonian for mesons. The spectrum of the renormalized 
Hamiltonian agrees with the one obtained with a smeared hyperfine interaction. The masses of the low-lying 
pseudo scalar and vector mesons are reasonably described within the model. 
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1. Introduction 

We address to the effective mass operator equa- 
tion of the tJ,-model for the qq Light-Front Fock 
state component of the meson with mass M: 



(1) 



virith the kernel 



U{k,p) = - 



4a 



2m^ 



Equation is mathematically not defined. 

The aim of this paper is to give Eq. a physi- 
cal meaning by renormalization, i.e., by applying 
to it the "subtraction method" , a renormalization 
scheme for nonrelativistic quantum mechanics of 
singular interactions developed earlier ||l| . 

This is an interesting problem because Eq.([T]) 
as proposed in Ref.||^, the tJ,-model, is an effec- 
tive Hamiltonian derived from Quantum Chro- 
modynamics, meant to describe the lowest Fock- 
state component of the Light-Front meson wave- 
function. It has been applied with reasonable 
success to the low-lying pseudo-scalar and vec- 
tor mesons, by using a different renormalization 
scheme namely by regularization and subsequent 
renormalization To the authors knowledge, 
fortunately nobody tried at that time in the past 
to solve the hyperfine spliting in the hydrogen 
atom using the Schroedinger equation directly. 



We are thus in the unique position to compare 
two drastically different schemes, both conceptu- 
ally and numerically, and verify that they agree. 
This strong statement stays at the very basis of 
renormalization ideas, that no matter the inter- 
mediate steps one performs to mathematically de- 
fine the initial equation (^, after renormalization 
all them produce the same physics. 

2. Notation 

For the purpose of presenting the subtraction 
method of Ref . , we introduce the notation be- 
low, and allow as well different quark masses and 
the relativistic phase space. 



(M^ + V + V^) 



\ip >-- 



(2) 



where the free mass operator of the quarks with 
masses mi and and m2 is Mq — Ei -\- E2. The 
individual energies are Ei — ^ mf + («=1,2) 
and k = \k\. The Coulomb- like effective potential 
is V and the hyperfine singular interaction is , 
which in the non-relativistic limit is the Dirac- 
delta. The matrix elements of these operators in 
momentum representation are given by: 

4m., ^, a 



<k\V\p>^~^ <k\x> -^<X\P>. (3) 



A 



< k\V \p> = < k\x> — < x\p> 



(4) 



where the total and reduced mass are denoted by 
nis — mi+m2 and m^ = mim2/ms, respectively. 
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A is the bare strength of the Dirac-delta hyperfine 
interaction. The mean four-momentum transfer 
is taken as = (fc — p)^. The phase-space di- 
mensionless function is 



1 



A{k) 



E1 + E2 
E1E2 



(5) 



The form-factor is x(fc) =< k\x >= A{k). 

3. Example: Dirac-Delta potential 

Many authors in the past have renormahzed 
the Schrocdinger equation with contact interac- 
tion (see e.g. |^). Here we want just to sup- 
ply the essence of the "subtraction method" of 
Ref.§. Let us solve Eq.(§) with only 



Uik,p) = X, 



(6) 



which, as said, makes Eq.(|l|) not well defined. 

The scattering amplitude for the potential (^) 
is the geometrical series 

t(M2) = [A-i -/(M2)]-\ 

solution of the scattering equation 

T(Af2) = V + VG''f^\K'P)T{M'^), 

for the scattering state of mass M . The Green's 
function of the free mass operator equation with 
outgoing wave boundary condition is 

The function 



dk 



1 



Af 2 _ 4^2 _ 4^2 _|_ 



diverges linearly! This is the mathematical prob- 
lem in Eq.([I]). 

How to give meaning to t(M^)? We use the 
renormalization idea. Suppose T(/i^) is known 
from experiment, then we rewrite t(M^) using 
this piece of data: 



{M^)=[T-\^,^) + I{^,^)-I{M^)]- 



(7) 



and know the subtraction of the divergence ap- 
pears! A closer look to 



dk 



(/i^ - 4m2 - 4A;2 -|- ie){M^ - Am? - Ak"^ + ie) 

shows that it is finite with /i being the subtrac- 
tion point! This is the essence of the "subtraction 
method" of Ref.ltl. 

The renormahzed form of the potential (||), V-^, 
has the bare strength written as a function of the 
renormahzed one XnifJ-^) = '''(m^) 



A 



1 



l + A7j(/.2)/(^2) 



(8) 



in which the physical input and the counter terms 
that subtract all the infinites in the scattering 
matrix at the mass scale /i are present [0. 

4. Renormalized Model 

The "subtraction method" exemplified in sec. 3 
is applied to the effective model defined by the 
mass operator of Eq.(j^). The scattering matrix 
comes from the solution of the scattering equation 
with the renormalized potential||l|,Q 



(9) 



where V-jz — V + V^. In finding the bound state, 
one could as well diagonalize the mass operator 

The renormalized Dirac-delta interaction is 
written formally as below 



- 



l + Ti{^?)G^^\^?) 



(10) 



where T^(/U.^) is the renormalized T-matrix of 
the Dirac-delta interaction, with matrix elements 
given by 

< p\Tii^l')\q>^ x{p)^ni^^')x{q)■ 

The solution of Eq. (|) is g: 

Tn{M^) = r^(M2) + |F > iTC(M') < F|, (11) 
where 

|F>= (l + r^(M2)G{,+)(M2)) \x>, 
<F\ =<x\ (Gi+\M')T''iM^) + l), 
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g{M) x\ [M' - Ml -V + ie] ' Ix >, 



and 



00 (M) xl [Af2 



Ml 



\x> 



The regular potential T-matrix, T^Af^), is the 
solution of the scattering equation (||) for V. 

The scattering equation with the renormalized 
interaction appears in a subtracted form [0, in 
which all the divergent momentum integrals are 
explicitly removed: 



Tn{M^)^T^^[^Ji'')[l 



(12) 



For a regular potential Eq.(12) is completely 
equivalent to the traditional Lippman-Schwinger 
scattering equation. 

We use the renormalization condition that at 
the pion mass, M = m^, the T-matrix, for mi = 
7712 — m„ = md has a bound-state pole, where 
t^{m'^) = 0. The choice /i = ni-^ implies that 

The physics described by the theory does 
not dependent on the arbitrary renormalization 
point, this imposes -^V^ — and qualifies the 
interaction as the fixed-point of this equation. 
The inhomogeneous term of Eq.jl^) runs as the 
subtraction point moves, according to the Callan- 
Symanzik equation 

The renormalized T-matrix (hi 



IS invariant un- 
der the change of fi to ^' , and thus -^t'ji{A'P) — 
0, and the renormalized strength runs as accord- 
ing to (m'2 ) = A-i {ii^)+giti) -g{p')-goit^) + 
goif^')- 

The bare strength is obtained by equating 
Eq.(||) to Eq.(p^), and using Atc(toJ), one finds 
Abore = rn'!^g{mTr)~^ , with the reduced mass 
mjy = 1/2 and ^/(m^) calculated for mi = m2 = 
rriu — m-^. With this, the pole of the T-matrix 



at the bound-state mass, Mb, is given by 
*niMi) = ^^("'-) - Gi^'^b) = , (14) 
for s-wave states with any quark mass. 



5. Comparing Renormalization Schemes 

Here we compare the results obtained with the 
Yukawa form of the smeared Dirac-delta inter- 
actions , Eq. ( |T5| ) , and the "subtracion method" , 
Eq.(|l^with A{k) = 1. Using Eq.(|l|), the pion 
mass, TOtt, of 140 MeV and the first excited s-wave 
state mass, m* , of 768 MeV were fitted with the 
parameters a = 0.763 and 77= 1148MeV see 
for TO = 406 MeV. The solution of Eq.(|l^), for 
TO^ = 140 MeV and a = 0.763 is to; = 766 MeV, 
in remarkable agreement with the previous result. 

Both renormalization methods have the same 
physical inputs. In one set of calculations, a was 
varied, with a fixed TO7r= 140 MeV. In the other 
set of calculations, nil= 768 MeV was kept fixed. 
In the model of Eq. (jlj) the value of 77 was fitted 
to the TOjr or m* for a given a. 

In figure 1, the results of to^ as a function 
of a for the two renormalization methods are 
shown. The agreement between the "subtraction 
method" and the smeared delta renormalization 
method is within few percent, which we relate 
to the rather drastically different methods. The 
values of 77 for a going to zero increase towards in- 
finite, to keep the ground state at the pion mass, 
while to; tends to the scattering threshold at 812 
MeV. For a increasing the values of rj decreases 
to keep TOtt fixed, and the to; which is Coulomb 
dominated, has to decrease, as we observe in fig- 
ure 1. The effect of the relativistic phase-space, 
Eq.(||), has been studied in Ref.([Q) and it is of 
the order of only few percents. 

The results for function of a for m* — 

768 MeV, are presented in figure 2. The thresh- 
old for zero pion mass occurs for a with the 
value about 0.75. The value of TOtt increases with 
a, corresponding to a decreasing binding energy, 
which means that the intensity of the short-range 
interaction, that dominates the ground state, 
diminuishes. In fact to keep constant to;, as 
the effective Coulomb interaction increases it de- 
mands a weaken short-range interaction. The cal- 
culation of TOtt with Eq.(|l5|) does not go beyond 
a — 0.97 because 77 vanishes and the mass of 768 
MeV of the excited state is reproduced with the 
effective Coulomb interaction. The "subtraction 
method" does not present the same limitation. 
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Figure 1. The mass TO*(MeV) is plotted versus 
a for a fixed = 140 MeV. The dashed curve 
gives results from Eq.(p^), the empty boxes from 
Eq.M) with A{k) = 1. 



Figure 2. The mass TO^(MeV) is plotted versus 
a for a fixed m* = 768 MeV. The solid curve 
gives results from Eq.(p^), the empty circles from 
Eq.(]l6|) with A{k) = 1. The two upper curves are 
the ones of fig. 1. 



6. Effective Meson Model 

Now, we use the "subtraction method" applied 
to the tJ,-model to calculate the mass gap be- 
tween the ground states of the pseudo-scalar and 
vector mesons, corresponding to (tt, p) [139,768], 
[494,892], (DO, [1865, 2007] and 
(B±,B*)[5279,5325] (experimental values of the 
meson masses within the square brackets in 
MeV). The ground state masses of the pseu- 
doscalar mesons comes from the solution of 
Eq.(|l^) with A{k) given by Eq.(|). The pion 
mass is fixed to its experimental value and the 
mass 7712 of one of the constituents quarks are 
varied with mi — 406 MeV and a constant. We 
choosed the values for a — 0.18, 0.4 and 0.5. 

The vector meson mass is associated to the sum 
TOi + TO2- In that sense, the vector meson mass 
does not have the contribution of the strong at- 
tractive hyperfine interaction, and the Coulomb 
effective attraction produces a binding energy too 
small compared to the mass, which we have de- 
sconsidered here. In figure 3, the results of the 
difference of the squared masses ml — m^^ of the 
ground state of the vector mesons and pseudo- 
scalar mesons as a function of the mass of the 
ground state of the pseudo-scalar mesons m^s are 
shown. The experimental results are reasonably 



described with a = 0.4. The value of a used to 
describe the data depends on mi. In our previ- 
ous model calculation, the value of 386 MeV was 
used and a was found to be 0.5 [^. The linear 
raising behaviour of the difference — m^^ with 
TOps observed in figure 3, is due to the saturation 
of the binding energy when m2 — > oo. It is not 
clear if a confining potential, not present in the 
model, can change such a trend. 

7. Conclusion 

The "subtraction method" [Q was applied to 
renormalize the ti-model||^ which contains an ef- 
fective Coulomb interaction and a hyperfine zero- 
ranged singular. We have compared with a dif- 
ferent renormalization scheme that make use of 
regularization and subsequent renormalization [^. 
The two drastically different schemes, both con- 
ceptually and numerically, agree. Here we provide 
one more simple example, that the physics of the 
renormalized theory does not recognize the in- 
termediate steps one performs to mathematically 
define the initial undefined theory. 

Acknowledgments: TF thanks to H. 
Leutwyler for suggesting the plot of figure 3. TF 
also thanks to H.C. Pauli for the warm hospitality 
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Figure 3. The difference ml — ■mp^{GeY'^) 
is plotted as a function of mps(GeV). Results 
from Eq.(p^), with a = 0.18 (dot-dashed curve), 
a = 0.4 (solid curve) and a — 0.5 (dashed curve). 
Experimental values (solid circles). 
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A. Equations in Momentum Space 

The effective model of Q corresponds to use 
the non-relativistic phase-space A(k) = 1 in 
Eq.(||) and a smeared delta-interaction of a 
Yukawa form: 



37r2 



2m^ 



^{p). (15) 



In the rcnormalizcd model for the Coulomb 
plus Dirac-delta interaction the bound state 
masses of the pion ground and excited states in 
s-wave, are found numerically from the zeroes of 



Eq.(14): 



dp 



47rp2 
lb) 



1 



1 



2 

r 

dq- 



q- 



dp 



^l-M§{p) 



VA{q) 
t^{p,q-mf) 



(m*2-A/2(p))(m*2_Af2(p)) 



{ml~M§ip))iml~M§iq)) 



0, (16) 



with m* the excited s-wave state mass. The 
free mass of the two quark system is Mo{k) = 
y/k'^ +mf + y/W+m^. The s-wave projected T- 
matrix of the Coulomb potential in Eq. (|l^) is 



i^(p,g;M^) 



dcos{e) <p\T^{M^)\q>; (17) 



which is the solution of 
4m a 



t^ip.q-M') 



1„ {p-iY 



Stt^pq ^A{p)A{q) 



8m I 
-—a / dp 
6n 



p' In 



p^A{p)A{p') - M2 (p) ' 



(18) 



the momentum space representation of s-wave 
projection of the scattering equation for the T- 
matrix T^(M2). 
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